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I. INTRODUCTION 

Logarithmic potential U{x) = g\og{x), freqnently appears as an effective potential in 
diffnsion models of biophysics and chemical physics. For instance, it is nsed to model 
interactions among polyelectrolytic polymers,^ it arises as an entropic term in free energy 
costs for nnzipping of DNA macromolecules,as an entropic potential stemming from 
interactions of colloids with walls of narrow channels,®^® and as a “potential” for momentum 
diffusion in optical lattices. 

Other prominent examples include long-range interacting particles,self-gravitating 
Brownian particles,diffusion of eigenvalues of random matrices known as Dyson’s Brow¬ 
nian motion,and interaction of tracers in driven lattice gasses.As for single-particle 
problems, diffusion in logarithmic potential has recently regained signihcant attention due 
to its intriguing slow relaxation^^^^'^’^^’^^ and rather nontrivial hrst-passage properties. 
Moreover, combining the logarithmic and the parabolic potentials^^’^®^^® one obtains exactly 
solvable models for stochastic resonance^® and for energetics of a Brownian particle in an 
asymmetric optical trap.®® 

In probability theory, the Brownian motion in logarithmic potential with a constant 
strength g is known as the Bessel process, which can be dehned as the distance to the 
origin of a ( 5 ^ -|- l)-dimensional Brownian motion.®^’®® In economy, scaling transformation 
of the Bessel process generate the so called Constant Elasticity of Variance model,®^’®® see 
Ref. 25 for more references in this direction. 

In the present work we consider the overdamped Brownian motion in the logarithmic 
potential with time-dependent strength 

U{x,t) = g{t)\og{x). (1) 

Then the position X(t) of the Brownian particle evolves according to the Langevin equation 

dX(t) = dW(t), (2) 

where D is the diffusion constant and W(t) denotes the standard Wiener process. Our 
primary interest is to characterize the first-passage time T, which is the time when the 
particle hits the coordinate x = 0 for the hrst time.®®’®^ The statistical properties of T 
are described by the survival probability S(t), also called persistence®®, which gives the 
probability that the particle has not hit the origin by the time t, S{t) = Prob{T > t}. 
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First-passage problems in chemical physics are nsed to model rates of diffnsion-limited 
react ions, exit kinetics of Brownian particles diffnsing out from biological micropores, 

or rates of any other process activated by diffusion. In our case one can imagine that Eq. (2) 
provides a rough model for the exit of a diffusing particle from a micropore with a time- 
dependent entropic barrier at the pore end (see Refs. 6-8, where it has been shown that 
the entropic potential is proportional to the logarithm of the pore width). At the same 
time, Eq. (2) can describe diffusion of a small mobile ion in the presence of a polyelectrolyte 
(long charged polymeric molecule).^ In the latter case, X(t) is understood as a distance 
of the ion from the charged polymer chain and the logarithmic potential can model the 
Coulomb potential energy of the small mobile ion [see Eqs. (4) and (15) in Ref. 1]. The 
strength g{t) of the logarithmic potential is then proportional to the density of charges on 
the polyelectrolyte. One can assume that this charge density depends on time due to the 
interaction of the polyelectrolyte with other small mobile ions. For instance the ions can 
“condense” on the polyelectrolyte reducing its charge density.^ 

To the best of our knowledge, there is no exact solution for the proposed model (2) 
when the potential strength g{t) depends arbitrarily on time. The present paper partially 
hlls this gap. We present exact asymptotic description of both the hrst-passage properties 
and the particle dynamics. Our main result is the exact asymptotic expression (12) for 
the survival probability S{t). This expression suggests that functions g{t) can be divided 
into three classes according to the resulting asymptotic decay of S{t) (see three dynamic 
regimes in Sec. IV B). The result (12) also determines the long-time asymptotic behavior 
of the probability density function (PDF) for the particle position, which yields asymptotic 
expressions for any moment (X"'(f)) as discussed in Sec. V. 

Finally it should be noted that, from a general perspective, the proposed model belongs to 
a class of hrst-passage problems with time-dependent conditions. Exact results in this held, 
even long-time asymptotic ones, are rather rare. Analytical progress is feasible for a free 
(no potential) Brownian motion in domains of varying size,^^^®^ for the Bessel process with 
drifting position of the singularity,for a Brownian motion in time-dependent piecewise 
linear potentials®^ (see also Ref. 59 for a related concept of “nonlinear relaxation time”), 
when the linear potential depends exponentially on time®'’^®^ and for a weak potential®® 
treated as perturbation. Our present hndings extend this list by a nontrivial example with 
a rather rich physical behavior. 
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We start the presentation by the precise dehnition of the model in terms of a Fokker- 
Planck eqnation supplemented by appropriate boundary conditions (Sec. II). Essential prop¬ 
erties of the Bessel process (constant g) are recapitulated in Sec. III. Sec. IV comprises 
our main results: exact asymptotic solution of the Fokker-Planck equation (3) (Sec. IV A) 
and classihcation of functions g{t) into three groups (three dynamic regimes in Sec. IV B). 
Asymptotic evolution of the particle position is discussed in Sec. V. Appendix describes the 
simulation algorithm used to check predictions of the asymptotic theory. 


II. DEFINITION OF THE MODEL 


Consider diffusion of a particle in a semi-inhnite one-dimensional interval (0, cxo) with an 
absorbing boundary at the origin. Initially, at f = 0, the particle is situated at xq, xq > 0. 
For f > 0 the particle diffuses in the logarithmic potential (1), hence its position evolves 
according to the Langevin equation (2). The boundary at the origin is perfectly absorbing, 
i.e., if the particle hits the origin, it is immediately absorbed with probability one. Thus the 
hitting time of the origin equals the time T of the hrst passage to the origin. 

The Fokker-Planck equation [corresponding to Eq. (2)] for the PDF of the particle position 
reads 


d d 


9{t) 

X 


p{x,t) 


( 3 ) 


The absorbing boundary at the origin is represented by the absorbing boundary condition 
p{0,t) = 0 and the initial condition is p(x, 0) = 6{x — xq). The survival probability S{t), 
that is the probability that T>t, follows from p(x, t) after the spatial integration 


S{t) = 


p{x, t)dx, 


( 4 ) 


Jo 

since S{t) equals the probability that, at time t, the particle can be found somewhere in the 
interval (0, oo). 


III. STATIC PHASE DIAGRAM {g = const.) 

Let us now recapitulate properties of the Bessel process (constant g), which will serve as 
a starting point in the following discussion. Derivations of discussed facts can be found e.g. 
in Ref. 23. 
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FIG. 1. Static phase diagram (constant g). In the region g/D < — 1 the logarithmic potential 
U{x) = 5 (log(x) is strongly repulsive as compared to the intensity of thermal fluctuations, therefore 
the particle never hits the origin. In all other cases {g/D > — 1), the survival probability decays 
to zero as t —)• oo. Its asymptotic decay is given in Eq. (7). Insets illustrate the potential for 
the individual regions of the phase diagram, U{x) is plotted for g = —1.5, g = —0.3 and g = I 
(assuming that D = 1). 


Dynamics of the model with a constant g exhibits two qualitatively different regimes. 
The master parameter which determines qualitative features of the dynamics is the ratio 
g/D. When g/D < —1, the potential is strongly repulsive as compared to the intensity 
of thermal fluctuations. Due to the strong repulsion, the Brownian particle never hits 
the origin x = 0 and its survival probability equals to one for any time. On the other 
hand, when g/D > —1, i.e., for the weakly repulsive potential [g/D G (—1,0)], for the free 
Brownian motion (^f = 0), and for the attractive potential {g > 0), the particle eventually 
hits the origin with probability one. Hence for g/D > —1 we have S'(f) —?■ 0 as t —?■ cxd. 
In Fig. 1, these dynamic regimes are summarized in what we call the static (constant g) 
phase diagram of the model. We can translate these conclusions for the ion diffusion in the 
presence of a polyelectrolytic molecule (see Introduction). When > 0, the small mobile 
ion and the polyelectrolyte possess opposite charges, thus it is not surprising that the ion 
will eventually be attracted to the polymeric chain. However, even when the ion and the 
polyelectrolyte bear charges of the same sign {g < 0), the ion can hit the polyelectrolyte 
provided the temperature is large enough and/or the charge density on the polyelectrolyte 
(which is proportional to g) is small, i.e., when 0> g/D> -1. 

In order to trace back the origin of such a nontrivial behavior notice that (for a constant g) 
the Langevin equation (2) will preserve its form after the “Brownian” rescaling of coordinate 
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and time, X(t) —)■ t Since the two competing tendencies - the diffusive 

spreading of the particle PDF and the drift of the particle away from the origin - scales 
with time in the very same way, it is the ratio of their amplitudes, g/D, which decides the 
type of resulting dynamics. 

In the present paper we will not discuss the case of the strong repulsion where gj D < — 1. 
We will focus on g such that g/D > —1, which includes both the region g/D > —1, where 
the survival probability decays to zero and the boundary case g = —D, where the survival 
probability equals to one for all times. The exact PDF of the particle position for g/D > — 1 
reads^^ 


PB(x,t) = 


Xq 

2Dt 


X 

Xo 


W-g/D) 


exp 


+ T (XXq\ 

ADt ) ’"\2Dt) 


(5) 


where the order of the modihed Bessel function of the hrst kind u depends on the potential 
strength, v = |(1 + g/D). The subscript “B” will mark all quantities related to the Bessel 
process. 

The survival probability S'B(t) is obtained from the PDF pB{x,t) after the spatial inte¬ 


gration, S'B(t) = pB{x,t)dx. For g/D > —1 the integral may be expressed in terms of 
the “upper” incomplete gamma function, F(a, 2 ;) = D~^e~*dt, as follows 


F 


SB(t) = 1 - 


Ul+W-D), 


■^0 

ADt 


( 6 ) 


r[i(i+ 9/-D)] 

For ADt ^ Xq, i.e., when the diffusion length is much larger than the initial distance of the 
particle from the origin, Eq. (6) simplihes to the power law 

.Xq/ _^_i (i+g/£i) 


SBit) 


r[i(3 + 9 /D)] 


t ^ OO. 


(7) 


Notice that the time-independent prefactor of the long-time asymptotic form (7) depends 
on the initial position xq. This long-time memory effect illustrates that the hrst-passage 
time T, which is characterized by SB^t), is a functional of the stochastic process X(t) (see 
reviews 36 and 65 for more details). Thus even in the long-time limit the asymptotic survival 
probability should bear information about the whole history of the process X(t). The latter 
fact is clearly visible in our asymptotic result (12), where the asymptotic survival probability 
depends on the integral of g{t). 

We remind that the PDF (5) is valid also for the boundary case which corresponds to 
g = —D. In this case, however, the PDF (5) is normalized to one for all times, i.e., S'B(t) = 1 
for all t. 
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IV. DYNAMIC PHASE DIAGRAM [g = g{t)\ 

If the strength of the logarithmic potential depends on time, g = g{t), then, to the best of 
our knowledge, the exact solution of the Fokker-Planck equation (3) is unknown. However, 
using the results for the Bessel process, it is possible to hnd the asymptotically exact solution 
of this equation for arbitrary g{t). In this section we present the theory and then discuss 
the asymptotic behavior of the survival probability S{t). The asymptotic expression for 
S(t) suggests that functions g(t) can be divided into three classes (three dynamic regimes 
in Sec. IV B) according to the resulting asymptotic decay of S{t). 


A. Asymptotically exact solution of the Fokker-Planck equation 


The central quantity for the following asymptotic theory is the PDF of the particle 
position given that the time of absorption T is greater than t, p{x,t\T > t). This PDF 
conditioned on nonabsorption before t, is dehned as 

p{x,t) 


p {x, t|T > t) = 


S{t) 


( 8 ) 


The approximation made below is motivated by a remarkable property of the conditional 
PDF of the Bessel process: its asymptotic behavior for any constant g is given by the 
Rayleigh distribution 

t ^ oo, (9) 


PB{x,t) X 


Ssit) 2Dt 

which can be verihed dividing the power series representation of PDF (5) by that of the 
survival probability (6). Thus regardless the actual strength of the logarithmic potential g, 
the conditional PDF is exactly the same as that for the free Brownian motion.®® This can be 
understood on physical grounds. When t is large, the PDF conditioned on non-absorption 
describes statistical properties of long-living trajectories. Such trajectories are typically 
situated far from the origin. At the same time, at large distances from the origin the force 
{—g/x) acting on the particle is vanishingly small and hence the long-time limit (9) does 
not dependent on the amplitude g. 

As a starting point of the approximative scheme we extrapolate the above observation 
to the case of time-dependent potential strength g{t). That is, in the long-time limit we 
approximate the conditional PDF p{x,t\T > t) by the Rayleigh distribution (9) and hence 
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for p(a;, t) we assnme the asymptotic approximation 

( 10 ) 


where S{t) is an unknown function yet to be determined. In order to derive S{t), which 
gives us the leading asymptotic behavior of the survival probability, we require that the 
Ansatz (10) for p(a:, t) satishes the exact Fokker-Planck equation (3) subject to the absorbing 
boundary condition p{0,t) = 0. 

If we insert the Ansatz (10) into the Fokker-Planck equation (3) it turns out that all 
x-dependent terms cancel. The hnal result is that the Ansatz (10) satishes the exact Fokker- 
Planck equation (subject to the absorbing boundary condition) if and only if the asymptotic 
survival probability S{t) satishes the hrst order ordinary diherential equation 


d ^ _2±ms^t). 

dt ^ ^ 2Dt ^ ^ 


( 11 ) 


Solution of Eq. (11) yields asymptotic behavior of the survival probability S(t), 

S{t) ^ A exp (^- ^ 2Df ^ 

where the lower limit to, to < t, is arbitrary since the constant (t-independent) term which 
originates from this limit can always be included into the t-independent prefactor A (the 
actual value of the prefactor cannot be determined from the above theory). In the whole 
paper we can simply set to = 0 except in Fig. 3, where to = 2 in order to avoid possible 
divergence caused by g(t), see the caption. 

Finally let us stress that the above derivation holds also for g{t) which lies on the phase 
boundary g(t) = —D, cf. Fig. 1. 


B. Three dynamic regimes 

The asymptotic survival probability (12) is the functional of the potential strength. As 
such it can exhibit a rather varied behavior for different functions g{t). Actually, Eq. (12) 
suggests that different g(t) can be grouped into three classes according to the type (regime) 
of the asymptotic decay of S{t). Specihcally, we distinguish 1. the regular regime, 2. the 
marginal regime, and 3. the regime of enhanced absorption. 

1. The regular regime is the only one where the asymptotic survival probability (12) can 
decay with time as the power law with a well dehned exponent. This regime occurs when 
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two necessary conditions are satisfied, (i) The first condition is that the potential strength 
tends to the limit g^o, —D < Qoo < oo, a.s t ^ oo. That is, we can express g{t) as 

g{t) = goc + e{t), (13) 

where the residne e{t) vanishes when t —)• oo. Then for S{t) we obtain 

S{t) ~ exp (^- , t^oo. (14) 

(11) The second condition is that the exponential in the above eqnation should converge to 

a constant as t ^ oo. To this end, the integral in the exponential must be bounded from 
above, [e{t')/t']dt' < oo, as t ^ oo. For example this is true for e{t) which decays to 
zero as an arbitrary power of time, e{t) ~ t~°‘, a > 0. Then the integral in the exponential 
tends to a constant as with increasing time and the exponential in Eq. (14) does 

not influence the leading asymptotic decay of S{t). Other prominent example is when e{t) 
oscillates around zero with a given period, say e{t) = Bsm{ut), for which [sin(ci;T)/T] dt' 
again converges to a constant. 

In the regular regime, when goo/D > —1, the survival probability decays as the power 
law jg^ exactly as that for the Bessel process (7) with g = g^o. On the 

other hand, even though the exponents in the both asymptotic formulas can be equal, the 
time-independent prefactors are in general different. The reason for this difference is that 
the constant A in Eq. (14) depends on the whole history of the process including the initial 
condition and the transient behavior of g{t). We are not able to determine the constant A 
within our asymptotic theory. In examples below, when we compare the analytical result 

(12) with simulated data, A is treated as the (only) fitting parameter. 

Besides the power law relaxation when goo/D > —1, the second possibility can occur 
within the regular regime. Namely, if g{t) tends to the phase boundary [goo = —D, cf. the 
static phase diagram in Fig. 1), the survival probability (14) converges to a nonzero value 
equal to A as t ^ oo. 

2. In the marginal regime the assumption (i) used to define the regular regime still holds 
and g{t) is again given by Eq. (13), whereas, the assumption (ii) concerning the fast decay 
of e{t) is no longer valid. Consequently, the exponential in Eq. (14) no longer tends to a 
constant as t —?■ oo. Instead, it behaves as a slow function of time, which modifies the 
asymptotic decay of S{t). A particularly simple example of the function e{t) which leads to 
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FIG. 2. Decay of the survival probability in the marginal regime for g{t) = —0.3 -|- B/log{t + 5), 
t € [0,10^], and for two values of B, B = —1 (the prefactor in the asymptotic survival probability 
(16) is A ~ 1.062) and B = -|-1 (^4 ~ 0.4021). The survival probability appears to be much larger 
for B = —1, as predicted by the logarithmic term in Eq. (16), which is shown in lower two panels. 
Other parameters used: xq = 1, D = 1; simulated data were averaged over 10® trajectories. 


the marginal behavior reads 


e{t) 


B 


t ^ oo. 


(15) 


log(i) ’ 

In this case, when Qoo/D > —1, the asymptotic survival probability tends to zero as the 
power law multiplied by the power of the logarithm; 


S{t) ~ At 2 ( 1 + 9 °°/'^) [log(t)] t 00 . 


(16) 


The logarithmic multiplicative term is clearly visible in simulated data, see two lower panels 
in Fig. 2. The sign of B signihcantly influences the absolute value of S{t) at hnite times. 
When B > 0, g{t) is always higher than its Qoo which leads to more rapid decrease of S{t). 
For B < 0, the value of g(t) at any hnite time is closer to the phase boundary as compared 
to goo- This results in much slower decrease of S{t) as demonstrated in the uppermost panel 
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FIG. 3. Slow logarithmic decay of S{t) according to Eq. (16) with g^o = —D, B = 1 and D = 1. 
In simulations (symbols) we have used g{t) = —1 -|- l/log(t), t € [2,10^], xq = 1; 10^ trajecto¬ 
ries were generated in total. The inset shows the difference between the simulated S{t) and the 
exact asymptotic behavior (16). Apparently the higher-order corrections to Eq. (16) also decay 
logarithmically with time when goo = —D. 

of Fig. 2. 

A particularly interesting situation arises when g(t) asymptotically approaches the phase 
boundary Qoo = —D. Then the power-law prefactor in Eq. (14) is absent and the asymptotic 
behavior of S{t) is determined solely by the decay of e{t). The resulting super slow decay 
of the survival probability for the logarithmic e{t) from Eq. (15) is demonstrated in Fig. 3. 
Eventually S{t) approaches zero, which is in sharp contrast to what we observe in the regular 
regime, where S{t) tends to a nonzero value equal to A, A < 1. 

3. The regime of enhanced absorption occurs when g{t) ^ oo as t ^ oo. In this case 
the potential becomes strongly attractive for large times (we move to the right in the static 
phase diagram in Fig. 1 as t —)■ cxo). Its strength can be decomposed as 

9{t) = go + r]{t), (17) 

where go, go > —D, denotes the value of g{t) at the initial time and ri{t) —)■ oo when t oo. 
For the asymptotic survival probability (12) we obtain 

S(t) ~ exp (^- j‘ j ^ oo. (18) 

In both the marginal and the regular regime S{t) decreases at most as the power of time 
[possibly multiplied by a slow function, see Eq. (16)] with the exponent being controlled 
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FIG. 4. Simulated S{t) (symbols) for g{t) = —D + t^, compared to the asymptotic result (20) (solid 
lines) for two different values of the exponent b, b = 0.5 (the fitted prefactor is ^4 ~ 1.141) and 
b = 1.5 {A Ri 1.015). Remaining parameters read D = 1, xq = 0.1. Simulated data were averaged 
over 10® trajectories. 


by the limit goo- Presently, since r]{t) grows with time, the overall asymptotic behavior is 
determined by the exponential in Eq. (18). The pre-exponential power law, which is now 
controlled by the initial valne plays a rather minor role. 

To illustrate the enhanced decay of -S'(t), let us assume a power-law increase of the 
potential strength 

9{t) = 90 + Bt’’, B >0, b> 0, (19) 

which yields the exponential asymptotic behavior 

S(t) ~ At--^^^+30/D) , t^oo. (20) 

\ 2Db J 

Value of the exponent b is decisive in the long-time limit, as demonstrated in Fig. 4, where 
the analytical result (20) is compared with simulated data for two values of b. On the other 
hand, the initial value go modihes just the pre-exponential factor. Its effect is however still 
well visible, see Fig. 5 which shows the results for different values of go- 

C. Failure of naive adiabatic approximation 

Assume that the potential strength g{t) varies rather slowly as compared to the diffusive 
spreading of the particle PDF. For such slow functions g{t), one can guess a straightforward 
approximation for S{t)- Namely, one can assume that at each time t the particle effectively 
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FIG. 5. Log-log plot of the pre-exponential factor in the asymptotic result (20) for three values of 
go- Simulation of S{t) (symbols) was performed for g{t) = gQ + Bt^ with B = 1 and b = 1/2. Other 
parameters used: D = I, xq = 0.5 when 5 o = 0 and go = D/2, and xq = 0.1 when g = —D. The 
fitted prefactor A reads A ^ 1.144 for go = —D, A « 0.390 for go = 0 and A « 0.186 for go = D/2. 
Simulated data were averaged over 10® trajectories. 


moves in a constant potential. This is the basic idea of the adiabatic approximation, which is 
freqnently a rather snccessfnl method how to obtain results for time-dependent problems.®^’®® 
In the spirit of the adiabatic approximation we can try to substitute the slow function g{t) 
for the constant g in the exact PDF for the Bessel process (5). The survival probability is 
then given by the spatial integral of such PDF. We thus obtain 


s^t) 


1 - 


F 


fliF+sWTFr 


( 21 ) 


which is the analogue of the survival probability (6) but now with time-dependent g{t). 
However, even for very slow functions g{t), this naive approximation gives different results 
as compared to the exact formula (12). Again the reason for this difference is that S{t) 
should reflect the complete history of the process, thus it is a functional of g{t) and not just 
a function of its instantaneous value. 
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V. MOMENTS OF THE PARTICLE PDF 


The asymptotic theory presented in Sec. IV A gives us also the leading asymptotic be¬ 
havior of all moments of the PDF p{x,t). From the dehnition of the conditional PDF (8) it 
follows that the nth moment of p{x, t) can be expressed as the product 


{X-{t))=S{t){X-{t))^^,, ( 22 ) 

of the survival probability S{t) and the nth conditional moment The latter is 

computed from PDF conditioned on nonabsorption before t, 

poo 

(X”(f))T>i= / x^p{x,t\T > t)dx. (23) 

Jo 

In the long-time limit the conditional moments are given by those of the Rayleigh distribu¬ 
tion, cf. the right-hand side of Eq. (9). Thus we have 

(X"(f))T>i ~ F(1 + n/2){4Dt)^/^, t ^ cx). (24) 


For instance for the hrst two moments of the particle position we obtain the asymptotic 
expressions 


(X(f)) ~ S{t)VM, (25) 

{X.\t)) ~ S{t)ADt, (26) 


where, for a specihc g(t), the asymptotic survival probability S(t) should be computed 
from Eq. (12) as discussed in the preceding section IV. Alternatively, one can obtain main 
asymptotic behavior of any moment (X"'(t)) from the moment-generating function 

Evaluation of the mean with respect to the PDF (10) yields 


~ >S'(f) 1 -|- ^\/ 7rDt e^^^^erfc ( —^\/Z^ 


(27) 


which, after expansion into a power series in gives the moments (24) multiplied by the 
survival probability S{t), cf. Eq. (22). 


VI. REGULARIZED LOGARITHMIC POTENTIAL 

Let us now consider a hrst passage time to the point a which is not exactly at the origin, 
0 < a < Xq. In order to reach a the Brownian particle has to overcome a finite (time- 
dependent) potential barrier, in contrast to the previous case a = 0, where the potential 
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FIG. 6. Log-log plot of Sa{t) for two values of g. The survival probability Sa{t) has been obtained 
by the numerical inversion of its Laplace transform (28) (symbols). Solid lines show only the 
first term of the asymptotic representation (29), dashed lines show both terms from (29). Other 
parameters used: D = 1, a = 0.5 xq = 0.6. 


barrier is inhnite. Such “regularization” of the logarithmic potential introduces an additional 
length scale which could signihcantly affect behavior of the survival probability. However, it 
turns out that the main asymptotic behavior of S{t) (now the probability that the particle 
has not reached the boundary at a; = a before f), as given by Eq. (12), remains valid also 
for the present regularized case. In the asymptotic formula (12) the regularization affects 
just the value of the constant prefactor A. 

New physical effects, which stem from the regularization of the potential, are not reflected 
in the time-dependence of the principal asymptotic part of S{t). Instead, an anomalous 
relaxation towards this asymptotic behavior is observed when a > 0 (as compared to the 
normal relaxation for a = 0). The relaxation is described by higher-order terms of the 
asymptotic expansion of *S'(f), which can be studied analytically for a constant g. In this 
case, the survival probability, say Sa{t), can be derived from the transition PDF (5) by a 
standard approach described e.g. in Refs. 69 and 70. Its Laplace transform, 


^a(s) 


1 ^ /xoy K^jxo^^s/D) 
s \aJ K^(av/s7D) 


(28) 


is valid for any g/D, i.e., for the whole phase diagram in Fig. 1. The order u of the modihed 
Bessel function of the second kind is determined by the ratio of the potential strength and 
the diffusion coefficient, u = ^{1 + g/D). In the region of the phase diagram which is 
close to the phase boundary, i.e., when g/D E (—1,1), the hrst two terms of the long-time 
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FIG. 7. The simulated survival probability for the regularized logarithmic potential (symbols) with 
strength g{t) = —D + t^ and for two values of b. Other parameters used: a = 0.3, xq = 0.4, D = 1. 
Simulated data were averaged over 10^ trajectories; asymptotic S(t) is given by (12). 


asymptotic series read^^'^^ 
xl'^ - 

r(i + 0 


Sa{t) 


h) +a 


21. - a 


2u 


r(i - 


4,Dt 


2u 


(29) 


^4:Dtj ' r(i — 2 i^) \r(i-i-z/) 

Contrary to this, when the potential is strongly attractive {g/D > 1) the second term of the 
asymptotic representation of Sa{t) wonld decay with time as 

First of all notice that indeed the leading asymptotic term in (29) is, up to the constant 
prefactor, identical to that given in Eq. (7) for the Bessel process. On the basis of the 
leading asymptotic term only, one could conclude that the regularized problem (a > 0) is 
rather similar to the case when a = 0. However, the second term in (29) reveals that for 
a hnite time t this conclusion is false. When g/D is close to —1 this second term vanishes 
extremely slowly with time. Thus for a successful quantitative comparison of Sa{t) (or of 
the simulated data) with the asymptotic representation (29) it is essential to include also 
this term (and possibly higher order terms). Fig. 6 demonstrates this fact for two values of 
g/D. Note that for g/D = —0.8, the second term is still visible at times t ~ 10^. 

The anomalous relaxation towards the main asymptotic behavior is absent when a = 0. 
In the limit a —)■ 0 all slow terms vanish and the asymptotic series of Sa{t) reduces to that 
for the Bessel process: 5*6(t) ~ cot~‘' (1 + ci/t + C 2 /t^ + ...), where the hrst-order correction 
vanishes (l/f)-times faster than the leading asymptotic behavior. 

Anomalous relaxation of the same kind is present when the potential strength depends on 
time. We have numerically simulated the process for several g{t) in each dynamical regime. 
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The simnlations snpport the above observation. Namely, when g{t)/D is close to the phase 
bonndary, we observed a slow relaxation similar to that in Fig. 6. When g{t)/D is far from 
the phase bonndary [e.g. g{t)/D > 1], then the asymptotic resnlt (12) shows satisfactory 
agreement with simulated data. Typical example of the latter situation is demonstrated in 
Fig. 7 (the regime of enhanced absorption), where g{t) departs from the phase boundary and 
grows as a power of time. Notice that for small times the simulated data differ signihcantly 
from the asymptotic curve. This is in sharp contrast to what we have observed for a = 0, 
cf. Fig. 4, where the asymptotic survival probability (12) agrees with simulation even for 
rather small times. 


VII. CONCLUDING REMARKS AND OPEN QUESTIONS 

Can a Brownian particle overcome the time-dependent logarithmic potential barrier (1) 
and hit the origin? If so, what are the statistical properties of the hitting (hrst-passage) 
time? When the potential strength is constant, the problem is exactly solvable and the origin 
can be reached only ii g/D > —1 a.s summarized in the static phase diagram in Fig. 1. In this 
case the survival probability decays as the power law (7). For a time-dependent potential 
strength the problem becomes more complex and, moreover, we have no exact solution of the 
Fokker-Planck equation (3) at our disposal. In the present paper we have presented simple 
long-time asymptotic theory which leads to the long-time solution (10) of this equation. 
The asymptotic PDF (10) factorizes into two terms. One term is the x-dependent Rayleigh 
distribution which is the conditional PDF for ensemble of long-living trajectories. Such 
trajectories are typically situated far from the origin, hence the conditional PDF contains 
no trace of the initial position xq and of the logarithmic singularity at x = 0. The other 
term, the asymptotic survival probability S{t) given by Eq. (12), bares information about 
both the initial position xq and the behavior of the function g{t') for all times t' before t. 
In Eq. (12) the initial condition is contained in the time-independent prefactor A [similarly 
as in Eq. (7) for a constant g]. The result (12) predicts rather rich asymptotic behavior 
of S{t) depending on concrete form of g{t). The three regimes of relaxation of S{t) can be 
identihed, see Sec. IV B. Having derived Eq. (12) one can return to the Ansatz (10) and 
thus obtain the asymptotic dynamics of the particle position, see Sec. V. Let us now turn 
to open questions which naturally stem from the presented results. 
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First of all a few further remarks on the natnre of approximation nsed are in order. Al¬ 
though the Ansatz for the long-time PDF (10) was motivated by the time-independent case, 
it turned out to be surprisingly successful when predicting the time-dependent behavior 
(see comparison with numerics for various g{t) in Figs. 2-5). This success stems from the 
fact that the differential equation (11) is closed, i.e., it contains no x-dependent terms, and 
thus the Ansatz (10) itself satishes the exact Fokker-Planck equation without any further 
approximation. On the other hand, on physical grounds one can expect that the key ap¬ 
proximation of the conditional PDF (8) by the Rayleigh distribution should at long times 
hold for a rather broad class of potentials. For instance for potential barriers localized near 
origin and for decreasing potentials U{x, t) = g{t)/x°', a > 0. The both examples yield force 
terms in the Langevin equation which decay faster than {—g/x) with increasing x. However, 
in these more general cases the resulting differential equation for the asymptotic S{t) is no 
longer closed and hence further approximations are needed. In future work it could be very 
interesting to extend the present theory to these potentials and/or to justify the present 
approximation rigorously, e.g., as an expansion in a small parameter. Such extensions of the 
present method could result in a unifying treatment of several rather nontrivial hrst-passage 
problems in time-dependent potentials so common in biophysics and chemical physics. 
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Appendix A: Simulation 

In order to verify asymptotic results we need the algorithm which fulhlls two requirements: 

(1) it treats exactly the neighborhood of the origin where the potential diverges; (ii) the 
average time step in the algorithm is large enough so that we could efficiently simulate the 
long-time regime. In order to meet both (i) and (ii), we have divided the code into two parts. 
Firstly, when the particle is far from the origin (when x > 1 in our code, this condition is 
satisfactory for all sets of parameters used in Figs. 2-5) we integrate the Langevin equation 

(2) by simple Euler-Maruyama method.Secondly, near origin (for x < 1 in our code) we 
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have generated a new position from the exact PDF (5) using the inverse transformation 
method/^ In the latter part the time-step At must be such that g{t) can be treated as 
constant in {t,t -|- At). For more details related to simulation of Bessel-like processes see 
Refs. 25 and 35. 
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